Duality in the Kondo model and perturbative approach to strong coupling theory 
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We investigate the duality in the Kondo model. Starting from the s-d model with the coupling 
constant J, the strong coupling model with the constant of 1/ J is identified. The model shows the 
unitary limit of the conductance, G = 2e 2 /h at zero temperature. The perturbation theory of the 
model gives qualitative agreement with the results of numerical renormalization group and Bethe 
ansatz near T = Tk- 
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Ever since Kondo elucidated the origin of logarithmic 
divergence of the resistivity in dilute alloy systems 0, 
various aspects of the Kondo effect have been investi- 
gated |2§, including recent experimental realizations in 

quantum dot systems HSIlEHIlH[l3- When 
the temperature or a characteristic energy scale is larger 
than the Kondo temperature Tk, the system is in the 
weak coupling regime while it is smaller, the system is in 
the strong coupling regime. In these regimes, the Kondo 
model exhibits duality. It was pointed out from the exact 
solution ^lj. In Refs. it was revealed for several 

quantum impurity problems, including S — 1/2 Kondo 
problem; In particular, the resistivity has calculated be- 
yond the Fermi liquid regime pl| , from the point of view 
of the duality. The duality also plays an important role 
in quantum dissipative systems 0, 0, 0], which are 
closely related with quantum impurity problems. It is 
possible to approach perturbatively for a strong coupling 
regime using a similar way in a corresponding weak cou- 
pling regime. 

We are presenting yet another duality for S = 1/2 
Kondo model with the coupling constant J. The moti- 
vation of this work comes from the following two points. 
The first one is a recent active interest on the Kondo 
screening cloud [l7L ll8L [TgL l2fj| . At low temperatures, 
conduction electrons and the impurity form the Kondo 
singlet state, accompanying a screening cloud. The na- 
ture of the cloud is still under debate, and accordingly 
many works have been investigated to detect the cloud 
in nano structures HH H3, H3, [2J, HE EE E3, IH, H^l - The 
second one is the notion of the strong coupling theory it- 
self. The Fermi liquid theories of the Kondo effect |30ll3l| 
have been developed and they successfully explain low 
temperature behavior of the Kondo effect Q|. Although 
it is stated that strong coupling theories describe the low 
temperature regime, simple questions have not yet an- 
swered clearly; What is the action (Hamiltonian) with 
the coupling constant of 1/J if it exists?, and how does 
the perturbation with respect to 1/ J work in the theory? 

We will explain that the concept of the screening cloud 
naturally leads to a duality in the Kondo problem. We 
will identify the model which has a very large cloud start- 
ing from the conventional weak coupling s-d model and 
show it has the coupling constant of 1/J. Electron trans- 
port of this model will be investigated, which shows the 



unitary limit of the conductance at zero temperature. 
A perturbative calculation with respect to 1/J with the 
poorman's scaling is performed, giving good qualitative 
agreement with the results of numerical renormalization 
group (NRG) [13 and Bethe ansatz [H around T = T K . 

Weak coupling model. — We start from the s-d model. 
The action S = So + Si is 
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where ipka is the conduction electron field with the mo- 
mentum k and spin a. We have introduced the surface 
field operator = Ylk^ka, which is the conduction 
electron field at the impurity site (x — 0). For simplicity, 
we have disregarded the potential scattering term. 

The two conduction electron fields, ip a and ipka can be 
treated independently if we introduce a Lagrange multi- 
plier A CT and add a constraint term to the action: 



J db^KQ) -Y^1>J\ +H, 



(4) 



The action can be written only by the surface field ip a 
after integrating over ip^ and X a : 



S 



E 



dwVv(w)G l {uj)ip (7 {oj) + Si, 



(5) 



where G(u) is the Green function of conduction elec- 
trons. At zero temperature, G(uj) — i/2 {2ttv) sgn(w) 
with the density of states v in the lead and the sign func- 
tion sgn(x). For finite temperature, it is convenient to 
introduce the Keldysh Green function (matrix) 



( 2mv[=£ + f(w)] 
[ ' ~ \ 2t^[/H - 1] 27m/- J 



2ixivf(iS) 

> /Ml 



(6) 



2 



where the matrix is defined on the Keldysh space [(1,1) 
component of the matrix gives the time-ordered Green 
function, for example.] 

Strong coupling model. — Now we show the duality of 
the Kondo model. Instead of integrating over A CT , we inte- 
grate over tp a and ipktr , and then the action is represented 
by A CT in the form: 



S 



duj\ rT (uj)g 1 (uj)\ <7 (u>) 



with 
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where we have used the identity of (<? ■ S)(a ■ S) = S 2 + 
io- S x S, and S x S = iS. The Green function g(oj) is the 
inverse of G(uj). For example, g(ui) — —2i/(2ni/) sgn(w) 
for zero temperature, and the Keldysh Green function is 



5M = 



fe[/M-i] 
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(9) 



Equation JJJ looks quite similar to the weak coupling 
model; It includes the s-d interaction as in Eq. J2J with 
the coupling constant proportional to 1/J instead of J. 
It also contain an additional "potential scattering" term 
which comes from the commutation relation among spin 
operators. We call this model as the strong coupling 
Kondo model, and define j = 16/(3J). 

Now we discuss the procedure has been made. The role 
of the Lagrange multiplier is to enforce the surface field 
■0,7 to locate at the impurity site; The interaction with 
the impurity is allowed only at the origin. As the tem- 
perature decreases, conduction electrons start to screen 
the impurity, accompanying the screening cloud, which 
is the sign of the Kondo effect. The size of the cloud is 
~ vf/Tk with the Fermi velocity vf of conduction elec- 
trons. It is estimated |24| to be the order of 1/jm in a 
recent experiment It is much larger than the size of 
the quantum dot. The large cloud means that conduc- 
tion electrons away from the origin can interact with the 
impurity strongly. Thus the surface field, which interacts 
with the localized spin, is no longer unnecessary to locate 
at the origin. As a result, ifj a — J2k fluctuates when 
the cloud is formed, and accordingly should be integrated 
out from the action. This means that the conjugate field, 
Ao- now describes the low energy physics instead. In this 
sense, the strong coupling model Q will describe the low 
temperature regime of the Kondo effect. In the following, 
this point is examined. 

Poorman's scaling. — Next we discuss the scaling of j 
of the strong coupling model using the poorman's scaling 
as in the weak coupling theory |34j . To begin with this, 
we includes the "potential scattering" term as the zeroth 



order action. Then the time-ordered Green function is 

1 (j/w 



1 + (j/2nu) 2 1 + 0727TI/) 



We perform the second order perturbation expansion 
with respect to j and then the RG equation reads 



dj 



dlogD 



(+4/(47T 2 i/)) 
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(ii) 



with an energy cut-off D. The second term in the time- 
ordered Green function does not induce the logarithmic 
contribution. 

There is an important difference in this RG; the sign of 
the right hand side is the opposite to the one in the weak 
coupling model. It comes from the difference in the sign 
of the time-ordered Green function. This result means 
that j becomes smaller as D — > 0; The s-d interaction is 
less prominent at low temperatures, and eventually van- 
ishes at T = 0. In terms of the weak coupling model, 
j = corresponds to J — ► oo, a strong coupling fixed 
point. Note that the scaling flow is formally the same as 
in the ferromagnetic Kondo model. The scaling flow al- 
lows to disregard j in the denominator in Eq. (|llfl which 
comes from the potential scattering term. Thus we dis- 
regard this term in the following. Then the solution is 



3(D) = 



Jo 



1 + 1/Wjo log d / D 1 



(12) 



where jo = j{do) with the initial cut-off do. 

To make a connection between the weak and strong 
coupling theories, jo and do should be represented by Jo 
and Do- The first condition is j(do)/v — 16 / '3(1 /uJ (do)) 
with vJ(do) = 1/ log(do/Tff). (do should be larger than 
Tk-) This results in 



3o(D)/v = 



1 



Wlo S (da/T K ) 



■ \og(d /D) 
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The value of do is still arbitrary, and the second condition 
is required. Equation H13(l diverges at D — di, which is 
a function of do- We impose the condition dd\/ddo = 
because the divergence should not depend on th e artifi- 
cial parameter do- This gives do = exp(y^3/167r)Tff ~ 
1.36Tr-. The two conditions eventually yield 



3(D)/u 



log(T*//J) 



(D < d ) 



(14) 



Since is 



with T£ = exp(27r- v /3/16)T i f ~ 15. 1T K 
larger than Tjc, j is a smooth function of D near T = Tk 
while J diverges at Tk- This provides a contrast between 
two perturbation theories as will be shown below. 

The expression of the current. — We shift our discus- 
sion to electron transport, and consider the model with 
the impurity attached to two leads a — L, R with a chem- 
ical potential difference \i. The Hamiltonian is Ho + H\ 
with H = J2 a .k,a^aka(^ ± A*/ 2 - Zktyako, [the sign 
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of +(-) is for L(R)], and Hi = JJ2 a , a > ■ aip a >. 
We have assumed fully symmetric constants between 
a = L, R for simplicity. The model is simplified, when 
we apply a unitary transformation between the con- 
duction electron fields, ip^ c = 1/ 'V2(ipL(k)<r + i>R(k)a) 

and ip( k)a = l/v / 2(V>L(fc) CT - i>R(k)a); The ip{k)a held 
is decoupled from the impurity, and it simply repre- 
sents free electron gas. A similar transformation is used 
in the Anderson model j^. We discuss the conduc- 
tance through the impurity. The current operator I(t) 
^[i>L{t)S • atpn{t) — ^)R{t)S ■ aipL(t)] is also rewritten by 



i = —j[ipS ■ dip — ips ■ dip]. 



(15) 



This term is added to the action in the form of 
Jdtrj(t)I(t). 

We repeat a similar procedure as above, rewriting the 
action using the Lagrange multiplier A Q . Applying the 
unitary transformation for X a , X v — 1/v / 2(Al + Xr) and 
A^, = 1/ v / 2(Al — \r), the two lead s-d model including 
the current generating term is given by 
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FIG. 1: The conductance as a function of \og 10 (T /Tk) cal- 
culated by various methods. The dots and thin solid line are 
for NRG [3^ and Bethe ansatz [3^ . respectively. The data 
is taken from Ref. |33|. The thick solid line is for the strong 
coupling perturbation theory. The broken line for the weak 
coupling perturbation theory |36|. 
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and qk = — i(2m/) tanh(w + /i/2)/2. In these matrices, 
we have represented them in the Keldysh rotation repre- 
sentation. 

The low energy action should be represented by X v and 
tp because the ip field participates in the screening while 
the ip field does not. Accordingly integrating over tp and 
A,/,, we obtain 



S = / du (X v ip) 



a b 
-b c 



16 f e 



dt r)(t)[X v ip - ipX v ] 
(17) 



with a = 



= A,b = 



i(fK-gK) 

2-KV 





The important 



result of this action is that the current generating term 
represented by X v and tp is independent on J and has no 
spin dependence although initially it depends on J and 
the spin. The lowest order of the conductance together 
with the result of the RG gives the unitary limit of G = 
2e 2 /h. 



Strong coupling perturbation theory. — We apply the 
perturbation theory with respect to 1/ J. The leading 
and next to leading order contributions give 



2e 2 



1 



3j 2 



\{2-Kvf 



(18) 



while G = = h -(2iruJ) ^ for the weak coupling 
model |3a| . We have disregarded the potential scattering 
term. To include the Kondo effect, the coupling constants 
should be replaced by the renormalized ones; J — > J(D) 
and j -> j(D). 

In Fig. 1, the conductance is plotted as a func- 
tion of log 10 (T/Tft:), calculated by various methods: the 
weak coupling perturbation theory j3|| , NRG j32| , Bethe 
ansatz |33j, and the strong coupling perturbation the- 
ory. The data of NRG and Bethe ansatz is taken from 
Ref. [3^. The perturbative strong coupling theory (the 
thick solid line) qualitatively coincides with the results of 
NRG and Bethe ansatz (the dots and thin solid line) near 
T < Tk- As T increases above Tk, the deviation from the 
results becomes larger. It is in contrast to the result of 
the weak coupling theory, where G deviates from the re- 
sults when T < IOTk- This difference describes the char- 
acter of the Kondo screening cloud. When T > IOTk, 
it is less prominent, and grows during Tk < T < IOTk- 
When T < T K , it is fully developed. 

Although the perturbation theory shows the unitary 
limit of the conductance, it does not reproduce the Fermi 
liquid picture at low temperatures; log(T) behavior still 
appears. This point is in contrast to the conventional no- 
tion of the strong coupling fixed point. The discrepancy 
comes from the limit of the perturbative RG approach. 
In the weak coupling theory, the perturbation with re- 
spect to J goes wrong as D decreases because it starts 
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from the doublet state rather than the singlet state. The 
same is true for the perturbation theory in the strong 
coupling regime. The Fermi liquid picture is captured by 
taking into account the spin singlet correlation appropri- 
ately from the beginning, for example applying a mean 
field theory for J 

Nevertheless, the strong coupling perturbation theory 
provides a simple approach to the crossover regime be- 
tween the Fermi liquid and the weak coupling theories, 
where has been treated only by more sophisticated or nu- 
merical methods. Furthermore, it elucidates that there 



is no transition around T — Tk, as has been shown. 

In conclusion, we have investigated the duality in the 
Kondo model to identify the strong coupling Kondo 
model with the coupling constant of 1/ J in a simple man- 
ner. The conductance calculated by the perturbation 
theory with the poorman's scaling qualitatively agrees 
with the exact results near T = Tk- 
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